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Abstract This paper has two parts: one about obser- 
vational constraints related to the empirical differential 
oxygen abundance distribution (EDOD), and the other 
about inhomogeneous models of chemical evolution, in 
particular the theoretical differential oxygen abundance 
distribution (TDOD). In the first part, the EDOD is de- 
duced from subsamples related to two different samples 
involving (i) N — 532 solar neighbourhood (SN) stars 
within the range, —1.5 <[Fe/H]< 0.5, for which the 
oxygen abundance has been determined both in pres- 
ence and in absence of the local thermodynamical equi- 
librium (LTE) approximation [Ramirez et al., Astron. 
Astrophys. 465, 271 (2007)]; and (ii) iV = 64 SN thick 
disk, SN thin disk, and bulge K-giant stars within the 
range, —1.7 <[Fe/H]< 0.5, for which the oxygen abun- 
dance has been determined [Melendez et al., Astron. As- 
trophys. 484, L21 (2008)]. A comparison is made with 
previous results implying use of [0/H]-[Fe/H] empiri- 
cal relations [Caimmi, Astron. Nachr. 322, 241 (2001b); 
Caimmi, New Astron. 12, 289 (2007)] related to (in) 372 
SN halo subdwarfs [Ryan and Norris, Astron. J. 101, 
1865 (1991)]; and (iv) 268 K-giant bulge stars [Sadler 
et al., Astron. J. 112, 171 (1996)]. The EDOD of the SN 
thick -|- thin disk is determined by weighting the mass, 
for assumed SN thick to thin disk mass ratio within the 
range, 0.1-0.9. 

In the second part, inhomogeneous models of chem- 
ical evolution for the SN thick disk, the SN thin disk, 
the SN thick + thin disk, the SN halo, and the bulge, 
are computed assuming the instantaneous recycling ap- 
proximation. The EDOD data are fitted, to an accept- 
able extent, by their TDOD counterparts with the ex- 
ception of the thin or thick -I- thin disk, where two 
additional restrictions are needed: (i) still undetected, 
low-oxygen abundance thin disk stars exist, and (ii) a 
single oxygen overabundant star is removed from a thin 



disk subsample. In any case, the (assumed power-law) 
stellar initial mass function (IMF) is universal but gas 
can be inhibited from, or enhanced in, forming stars at 
different rates with respect to a selected reference case. 
Models involving a strictly universal IMF (i.e. gas nei- 
ther inhibited from, nor enhanced in, forming stars with 
respect to a selected reference case) can also reproduce 
the data to an acceptable extent. 

Our main conclusions are (1) different models are 
necessary to fit the (incomplete) halo sample, which 
is consistent with the idea of two distinct halo com- 
ponents: an inner, flattened halo in slow prograde ro- 
tation, and an outer, spherical halo in net retrograde 
rotation [Carollo et al.. Nature 450, 1020 (2007)]; (2) 
the oxygen enrichment within the inner SN halo, the SN 
thick disk, and the bulge, was similar and coeval within 
the same metallicity range, as inferred from observa- 
tions [Prochaska et al., Astron. J. 120, 2513 (2000)]; 
(3) the fit to thin disk data implies an oxygen abun- 
dance range similar to its thick disk counterpart, with 
the extension of conclusion (2) to the thin disk, and the 
evolution of the thick -I- thin disk as a whole [Haywood, 
Mon. Not. R. Astron. Soc. 388, 1175 (2008)] cannot 
be excluded; (4) leaving outside the outer halo, a fit to 
the data related to different environments is provided 
by models with a strictly universal IMF but different 
probabilities of a region being active, which implies dif- 
ferent global efhciencies of the star formation rate; (5) 
a special case of stellar migration across the disk can 
be described by models with enhanced star formation, 
where a fraction of currently observed SN stars were 
born in situ and a comparable fraction is due to the 
net effect of stellar migration, according to recent re- 
sults based on high-resolution iV-body + smooth par- 
ticle hydrodynamics simulations [Roskar et al., Astro- 
phys. J. Lett. 684, L79 (2008)]. 
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1 Introduction 

The synthesis of a-elements (such as O, Ne, Mg, Si, 
S, Ca) takes place within SnII progenitors, via reac- 
tions of the type, E + a ^ E' + 7, where E, E', are 
generic elements (E lighter than E'), a is a helium nu- 
cleus, and 7 a photon. Accordingly, a-elements are 
primary elements, for which the yield is not sensitive 
to the initial abundance of carbon or other heavy ele- 
ments in the progenitor, because they are synthesized 
directly from hydrogen and helium. The abundances of 
a-elements as a function of the metal content provide 
crucial information about the stellar initial mass func- 



different dynamical signature is exhibited by the Galac- 
tic spheroid and disk, which makes evidence for a dis- 
tinct halo-bulgc and thick disk - thin disk collapse, 
as shown by different em pirical distributions of spe - 
cific angular momentum (|Wvse and Gilmord (Il992h : 
Ibata and Gilmord (|l995l )V 



tion ( I MF) a n d formation histo r y (e.g . , iTinslevI (119801): 
Pagell (Il989f): IProchaska et al.1 (|2000l ): iRamirez et al ' 



Concerning metal abundance distributions in SN 
disk stars, it has been recognized that earlier re- 
sults based on spectral type selection, were biased 
against stars with solar metallicity or higher, with re- 
spect to recent results based on colour selection (Hay- 
wood (2001, 2006). Oxygen abundance distributions 
in SN disk stars, or thick disk and thin disk sepa- 
rately , can be deduced fro m recent direct determina- 



( 20071 ): iMelendez et al.l (|2008l )). and make a main con- 
straint on models of chemical evolution. Among a- 
elements, oxygen has a special role, in that it is the 
most abundant element in the universe, after hydrogen 
and helium. 

Unfortunately, oxygen is more difficult than iron 
(which is not an a-element) to detect, and empir- 
ical relations were used to e xpress the form er as a 



tions (jRamirez et al.l (|2007f ) ) . The same cannot be done 
for halo and bu lge stars, with t h e exce pt ion of restricted 
sampl es (e.g., Ramirez et al.l ( 2007t ): iMelendez et al 



functi o n of the latter ( e.g., iBarbuv ( 1988 ): 



Carretta et al 



Page! 



(2000) 



( 2008f l). and some empirical [0/H]-[Fe/H] relation must 
be used (e.g., Caimmi ( 2007^ . hereafter quoted as C07). 
Recent analysis on high-resolution infrared spectra of 
giant stars, show no chemic al distinction between the 
bulge and the SN thick disk |Melendez et all |2008l) ). 

Inhomogeneous (i.e. implying inhomogeneous star 
formation) models of chemical evolution succeed in 
both providing a solution to the G-dwarf problem and 



( 1989l):lAbia and Rebolold 19891): 

Gratton et al.l (j2000l ): Takad a-Hidai et al.' ('2001)I llsraelian efepalc^ ucing substantial scatter exhibite d by the em- 
( 2001 ): for a review refer to|Barbuy ct al. (200l|). Only pirica l S N age-m e tallici ty relation ( e.g., Malinie et al 



recently, direct oxygen abundance determinations on 
solar neighbourhood (SN) thick disk, thin disk, and 
halo star jjj, and bulg e stars have been perfor med (e.g., 
Ramirez et al] (|2007l) : IMelendez etall (|2008l )). Differ- 
ent oxygen abundance at any given iron abundance 
range in common between two selected environments 
implies related stars were formed from different mix- 
tures of gas, presumably at different epochs and/or 
locations. Oxygen abundance patterns and their sys- 
tematic diversities between environment stars thus con- 
strain important information about the assembling and 
evolution of the Galaxy. 

In particular, the chemical abundance of thick disk 
stars suggests a similar history to those of metal- 
rich ([Fe/II]« —1.3) halo stars, and the thick disk 
abundance patterns show excellent agreement with the 
chemical abundances observed in metal-poor ([Fe/H]^ 
—0.4) bulge stars, suggesting the two populations were 
formed from the same reservoir at a common epoch, re- 
gardless of the exact physical proc ess involved in their 
formation ( Prochaska et al.l |2000l )). Additional sup- 
port to the last concl usion is provided by recent find- 
ings ( Melendez et al ] (jioos)). On the other hand, a 



( 1993h : ICaimmil (|2000h . ICaimmil ( 2001a ) 
quoted together as COO: IC"''^""'* 'It^'^itJ^ 



hereafter 



2001br ). hereafter 



^aimmi 

quoted as COl; C07: iCaimmi (200^, hereafter quoted 
as C08). The current paper aims to investigate if in- 
homogeneous simple models of chemical evolution with 
strictly universal IMF (implying gas is neither inhib- 
ited from, nor enhanced in, forming stars with regard 
to different environments) provide an acceptable fit to 
the oxygen abundance distribution in the SN thick disk, 
SN thin disk, SN halo, and bulge, for which new data 
are available, even if biased towards low metallicities 
or incomplete. Though the existence of a universal 



IMF ha s recently been questioned (e.g., I Bailer o et al 



( 2007al) ). it appears as a natural consequence of simi- 
lar (but not necessarily shared) chemical evolution his- 
tories experienced by the m etal-rich SN halo and t he 

(|200nh ). 



SN thick disk on one hand ( Prochaska et al 



and the bulge and the SN thick disk on the other hand 
( Melendez et"al1 |2008l )). 

The oxygen abundance is deduced from a sample 
(N = 523) of SN disk and SN halo stars jRamirez et al 



^In general, SN stars are to be intended as "passing through the 
solar neighbourhood" (e.g.. lRvan and NorrisI lll99ll) ). 



( 20071) ). lereafter quoted as the Ra07 sample, from 
which four subsamples are extracted following the pre- 
scriptions of the authors, related to thick disk {N = 
133), thin disk {N ^ 310), halo {N = 28), and im- 
certain population {N — 51) stars, hereafter quoted 
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as the RaK07, RaN07, RaH07, and RalJOT subsample, 
respectively. The RaK07 subsample is biased against 
low ([Fe/H]~ —1) m etallicities, and the R aH07 sub- 
sample is incomplete ( Ramirez et al.l ( 2007 )). For the 
bulge and the halo earlier results (C07) are also used, 
where oxygen abundance is deduced from iron abun- 
dance using empirical relations, with regard to a sam- 
ple ( N = 268) of bul ge K giants in Baade's win- 
dow dSadler et all (|l996h ). hereafter quoted as the Sa96 
sample, and a sample (N = 372) of kinematic ally se- 
lected halo subdwarfs ( Ryan and Norrij ( 199l| )), here- 
after quoted as the RN91 sample. 

The oxygen abundance is also deduced from a more 
restricted sa mple (N = 68) of SN disk, SN halo, and 
bulge stars ( Melendez et al. ( 2008h ). hereafter quoted 
as the Ma08 sample, which is subdivided into four sub- 
samples, related to the thick disk {N — 21), thin disk 
{N = 24), halo (iV = 4), and bulge {N = 19), hereafter 
quoted as the MaK08, MaN08, MaHOS, and MaB08 
subsample, respectively. All objects have similar stellar 
parameters but cover a broad range in metallicity. Oxy- 
gen and iron abundances were determined from a stan- 
dard one-dimension (ID) local thermodynamic equilib- 
rium analysis. Systematic errors were minimized using 
a homogeneous and differential analysis of the bulge, 
halo, thick disk, and thin disk stars. For further details 
refer to the parent paper ( Melendez et al. ( 20081 )). 

The value of solar oxygen abundance also affects 
the [0/Fe]-[Fe/H] relation. Recent determinations 



is determined in Section [2] for RaK07, RaN07, RaH07, 
RaU07, and MaK08, MaN08, MaB08, subsamples. In 
addition, a putative EDOD is derived for the SN thick 
-f thin disk by weighting the mass of the related subsys- 
tems, in Section [31 A comparison with the predictions 
of simple inhomogeneous models, is given in Sectional 
The discussion and related implications for the forma- 
tion of the Galaxy are the subject of Section [S] Some 
concluding remarks are reported in Section [6l 



2 The data 

2.1 The empirical differential oxygen abundance 
distribution 

While observations are related to logarithmic num- 
ber abundances, [E/H]=log(E/H) — log(Eo/HQ), mod- 
els of chemical evolution deal with mass abundances, 
Ze = M^/M, where E is a generic element heavier 
than He, Me and M are the total mass under the form 
of the element, E, and of a ll elements, respective l y. Th e 
following relation dPagell (Il989l) : iMalinie et all |l993l) ; 
Rocha-Pinto and Maciell (|l996l) : COO, COl, C07, C08): 



log (j> = log 



Zq 
{Zo)q 



(1) 



point t owards lower val u es (e. g .. lAllende-Prieto et 8,1 



( 20011): ISofia and Meved (|200lh : lAsplund et al.1 (|2004l) : 
Melendea (120041) 1. b ut the q uestion is still under debate 
(e.g., Landi et al. ( 2007); Socas-Navarro and NortonI 



( 2007 )): ICenteno and Socas-Navarrol(|2008l ): ICaffau eta H 
( 20081 )). A solar ox ygen mass abundance, (.^n )(Ti = 
0.0056, deduced from lAllende-Prieto et all (|200lh . shall 
be used to preserve comparison with earlier results 
(COl, C07, C08). 

Due to the above mentioned uncertainties, it has 
been preferred to deal with less recent iron abun- 
dance determinations to preserve comparison with ear- 
lier work (COl, C07, C08), where direct oxygen abun- 
dance determinations are still lacking. More recent, 
statistically more significant and less well scrutinised 
de terminations e.g., th e Geneva-Copenhagen survey 



(|Nordstr5m et all ( 20041)) for the SN dis k and the Ham- 
burg/ESO survey (|Schorck et all (|2008l )l for the metal- 



poor SN halo, are related to iron abundance instead 
of oxygen abundance as in e.g., Ra07 and Ma08 sam- 
ples, and for this reason shall not be used in the current 
attempt. 

The empirical, differential, oxygen abundance dis- 
tribution (EDOD), V = A7V/ (7VA(/i), log0 =[0/H], 



holds to a good extent. For a formal derivation refer to 
earlier work (C07, Appendix A)^ 

Fol lo wing previous at tempts (Pagel (^19891) : lMalinie et al.t 
( I993I ): Rocha-Pinto and Maciel (1996. )), the compari- 
son between model predictions and observations shall 
be performed using the differential instead of the cu- 
mulative m etallici t y dist ribution, as it is a more sen- 
sitive test ( Pagel ( 19891 )) an d allows direct compari- 
son be tween different samples (jRocha-Pinto and Maciel 
( 1996[ )). The occurrence of a sensitivity error in per- 
forming observations, precludes the use of a proper dif- 
ferential notation in dealing with the EDOD. This is 
the reason why the bin length cannot be lower than the 
sensitivity error, and differential ratios i.e. first deriva- 
tives, dN / d(/), must be replaced by increment ratios, 
AN/A(f>, where is the bin length, AiV the number 
of sample objects within the selected bin, and N is the 
total number of sample objects. 

Accordingly, the EDOD in a selected class of objects 
is defined as: 



V'(0T A^^ 



log 



NAcj) 



expio [0/H] 



(2a) 

(2b) 
(2c) 
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where in general, exp^ defines the power of basis, ^ 
(in particular, exp defines the power of basis, e, ac- 
cording to the standard notation), and the incremen t 
ratio, AN/ Acj ) , used in earlier attempts (jPagell (jl989l ): 
Malinie et al. ( 1993[ l) has been replaced by its normal- 
ized counterpa rt, AN/ (NAcj)), used in more recent 
investigations ( Rocha-Pinto and Maciell (1996); COO, 
COl, C07, COS), to allow comparison between different 
samples. The uncertaint y on AA^" has been evalu ated 
from Poisson errors (e.g., Rvan and NorrisI ( 199ll )). as 
CAJV = {ANy^^, and the related uncertainty in the 
EDOD is (e.g., COl, C07, COS): 



A'^lp = 1^/; - -0^1 ^ log 



It 



(AiV) 



1/2 



AN 



AN T {AN) 
= log ^ 



1/2 



(3a) 
(3b) 



where "0" — > —oo in the limit AA'^ — > 1. For further 
details refer to earlier work (COl). 

In the following tables and figures, the bin sizes in 
normalized oxygen abundance, 0, shall correspond to 
uniform bin sizes in [0/H] (or [Fe/H] if an empirical 
[0/H]-[Fe/H] relation is used), which implies non uni- 
form bin sizes, Acj) = 2A'^(f), as shown in Table [TJ 

2.2 Oxygen abundance distribution in subsamples 

The Ra07 sample is considered, where oxygen abun- 
dance is directly inferred both in presence and in ab- 
sence of the lo cal thermodyna r nic eq uilibrium (LTE) 



approximation (jRamirez et al.l (j2007l) . Table 6, avail- 



able at the CDS). The LTE approximation well holds 
deep in the stellar photosphere, and a continuous tran- 
sition occurs up to complete non-equilibrum (NLTE) 
high in the atmosphere. For further details refer to 
specialized textbooks (e.g., Gray ( 2005h . Chapter 6). 

The metallicity range is —1.5 < [Fe/H] < 0.5, 
but the sample is biased against lower metallicities, 
[Fe/H]'^' —1. On the other hand, it is expected to 
be unbiased at higher metallicities, for the following 
reasons: (i) more than 90% of the sample stars are 
located within 100 pc from the sun, and (ii) the re- 
lated differential iron abundance distribution looks sim- 
ilar to it s counterpart deduced fro m different SN dis k 
samples (Wvse and Gilmorel (|l995h: I.Torgensenl (120001 ): 



HavwoodI (|200lD : 



Fuhrmannf ( 200Sl )). and the same 



holds for SN thick disk and thin disk subsamples, when 
comparison may be performed. The EDOD has been 
deduced from the Ra07 sample in an earlier attempt 
(COS), both in presence and in absence of LTE approx- 
imation. 

Under the assumption that sample objects belong to 
three stellar populations (thick disk, thin disk, halo). 



each with a Gaussian velocity distribution, the proba- 
bility of a star being located in one of the components, 
may be expressed as a function of the Galactic space 
velocities. Pi — Pi{U,V,W), where i = 1,2,3, for the 
thin disk, the thick disk, and the halo, respectively; a 
Pi > 0.7 constraint is set to determine t he membership 
of sam ple objects to the i population ([Ramirez et al 
(l2007h V 



Using the above mentioned prescriptions, the Ra07 
sample has been subdivided into four subsamples: 
RaN07, RaK07, RaH07, RaU07, where Pi > 0.7, 
P2 > 0.7, P3 > 0.7, P^ < 0.7, i = 1,2,3, respec- 
tively. The RaH07 subsample is by far incomplete, 
as [Fe/H]> —1.5 for sample objects and any attempt to 
trace a halo abundance trend using only the metal-rich 
tail is dubious (Ramirez et al., 2007). With this caveat 
in mind, the RaH07 subsample shall be considered to- 
gether with the RaU07 subsample. 



On the other hand, it can be seen ([Ramirez et al 



( 2007t l) that a similar [0/Fe]-[Fe/H] empirical relation 
is exhibited by RaK07 and RaH07 subsamples: 



[0/Fe] =ao + 6o[Fe/H] 



(4) 



where ao = 0.370 T 0.027, 60 = -0.121 =F 0.043, Nq = 
43, (TO = 0.065, for the RaK07 subsample; ao = 0.3SST 
0.049, 60 = -0.048 T 0.071, Nq = 12, ao = 0.072, for 
the RaH07 subsample; in addition, A^o is the number of 
stars used for the fit and ao the random scatter around 
the mean fit. 

The MaOS sample is also considered, where oxygen 
abundance is deduced from a standard ID local ther- 
modynamic equilibrium analysis, and systematic er- 
rors are minimized using a homogeneous and differen- 
tial anal ysis of s t ars be longing to different populations 
(Melende zet al.l |2008h ). The MaOS sample is made 



of four subsamples: MaKOS, MaNOS, MaHOS, MaBOS, 
within the metallicity range, —1.7 < [Fe/H] < 0.5. 

A similar fit, expressed by Eq. ((4|), is shown using 
MaBOS and MaKOS subsamples, where ao = 0.41, 
bo = -0.02, No = 19, ao = 0.09, for the MaBOS 
subsample, and ao = 0.39, bo = — C j.Ol, A^o = 21 
ao = 0.09, for the MaKOS subsample ([Melendez et al. 
( 200S[ )). It can also be seen that oxygen abundance 



patterns exhibited by the MaHOS subsample are very 
close to their count erparts related to the MaBOS and 
MaKOS subsamples ([Melendez et al.l ([200Sh . Fig. 2). 

In conclusion, a similar trend for the [0/Fe]-[Fe/H] 
relation with respect to the SN thick disk, the bulge, 
and the SN halo, is deduced from both Ra07 and MaOS 
samples, in the metallicity range where the thick disk 
is unambiguously defined. 
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2.3 Oxygen abundance distribution in SN thick disk 



Both RaK07 and MaK08 subsamples are considered to 
infer separate EDODs related to the SN thick disk, 
for comparison with model predictions. The oxygen 
abundance range is -0.60 < [0/H] < 0.45 (LTE), 
-0.80 < [0/H] < 0.45 (NLTE), for the RaK07 sub- 
sample, and -1.40 < [0/H] < 0.30 for the MaKOS 
subsample. 

The EDOD related to the RaK07 subsample is listed 
in Table [H and plotted in Fig.[Il upper (LTE) and 
lower (NLTE) left, where the dashed vertical bands 
correspond to [Fe/H]= —1 and related uncertainty, 
deduced from Eg. (HI) in dealing with the thick disk 
([Ramirez et al.l (|2007t )). 

The EDOD related to the MaKOS subsample is hsted 
in Table [3] (left) and plotted in Fig.[2] (upper left). 

It can be seen that EDODs deduced from RaK07 and 
MaKOS subsamples show qualitative agreement, leav- 
ing aside the low-metallicity ([Fe/H]< —1) tail where, 
on the other hand, oxygen abundance bins are poorly 
populated and the statistical relevance is weak. 



Table 1 Logarithmic and numerical bins in oxygen abun- 
dance. Upper and lower values for each bin are denoted as 
[O/H]^ = log<^^, according to Eqs. ([TJ and (Hb)). Logarith- 
mic bins which differ by 1 dex are related to numerical bins 
which differ by a factor 10, which makes an (up and down) 
endless recursion of the table 



[0/H]- [0/H]^ 






05 





15 


1 


2673 


EH 


rO 


1 


4526 


E- 


1 





15 





25 


1 


5954 


EH 


-0 


1 


S2S7 


E- 


1 





25 





35 


2 


00S5 


EH 


-0 


2 


3022 


E- 


1 





35 





45 


2 


52S6 


EH 


hO 


2 


S9S3 


E- 


1 





45 





55 


3 


1S33 


EH 


hO 


3 


64SS 


E- 


1 





55 





65 


4 


0075 


EH 


-0 


4 


5935 


E- 


1 





65 





75 


5 


0451 


EH 


hO 


5 


7S29 


E- 


1 





75 





S5 


6 


3514 


EH 


hO 


7 


2S02 


E- 


1 





S5 





95 


9 


9960 


EH 


-0 


9 


1653 


E- 


1 





95 


1 


05 


1 


0066 


EH 


hi 


1 


153S 


E- 





1 


05 


1 


15 


1 


2673 


EH 


^1 


1 


4526 


E- 






2.4 Oxygen abundance distribution in SN thin disk 

Both RaN07 and MaNOS subsamples are considered 
to infer separate EDODs related to the SN thin disk, 
for comparison with model predictions. The oxygen 
abundance range is -0.35 < [0/H] < 0.60 (LTE), 
-0.45 < [0/H] < 0.45 (NLTE), for the RaN07 sub- 
sample, and -1.00 < [0/H] < 0.30 for the MaNOS sub- 
sample. 

The EDOD related to the RaN07 subsample is hsted 
in Table[l]and plotted in Fig.[Tl upper (LTE) and lower 
(NLTE) right, where the dashed vertical bands cor- 
respond to [Fe/H]= —1 and related uncertainty, de- 
duced from Eg . (IH) in dealing with the SN thin disk 
( Ramirez et"alj (|2007ll ) . 



The EDOD related to the MaNOS subsample is listed 
in Table [3] (right) and plotted in Fig. [2] (upper right). 

It can be seen that EDODs deduced from RaN07 and 
MaNOS subsamples show qualitative agreement, leaving 
aside the low-metallicity ([Fe/H]< —1) tail where, on 
the other hand, oxygen abundance bins are poorly pop- 
ulated and the statistical relevance is weak. 

2.5 Oxygen abundance distribution in the bulge 

The MaBOS subsample is considered to infer the EDOD 
related to the bulge, for comparison with model predic- 
tions. The oxygen abundance range is —1.3 < [0/H] < 
0.5. 

The EDOD related to the MaBOS subsample is hsted 
in Table [5] and plotted in Fig. [5] (lower right) . 



Table 5 The empirical, differential oxygen abundance dis- 
tribution (EDOD) in the bulge, deduced from the MaBOS 
(A'^ = 19) subsample 



1 BULGE 
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2 
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-4 
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Table 2 The empirical, differential oxygen abundance distribution (EDOD) in the solar neighbourhood (SN) thick disk, 
deduced from the RaK07 subsample {N = 133) both in presence (LTE) and in absence (NLTE) of the local thermodynamic 
equilibrium approximation 
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Table 3 The empirical, differential oxygen abundance distribution (EDOD) in the solar neighbourhood (SN) thick disk 
(left) and SN thin disk (right), deduced from the MaN08 {N = 24) and the MaK08 {N = 21) subsample, respectively 
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Fig. 1 The empirical, differential oxygen abundance distribution (EDOD) related to the RaK07 (left panels) and RaN07 
(right panels) subsample, both in presence (LTE, top panels) and in absence (NLTE, bottom panels) of the local ther- 
modynamic equilibrium approximation. The theoretical, differential oxygen abundance distribution (TDOD) related to 
inhomogeneous models of chemical evolution for the solar neighbourhood (SN) thick disk (case DB15, left panels) and the 
SN thin disk (case DN, right panels) are plotted as full curves, and as dotted curves in the remaining panels. Models are 
fitted to NLTE data. The dashed vertical bands correspond to [Fe/H[]= —1 and related unce rtainties, deduced from Eq. ^ 
in dealing with the thick (left panels) and the thin (right panels) disk (jRamirez et al.l lj2007h 



It can be seen that EDODs deduced from the MaN08 
subsample and the Sa96 sample, where [0/H]-[Fe/H] 
empirical relations have been used (C07), show quali- 
tative agreement, in spite of the poorly populated oxy- 
gen abundance bins, for which the statistical relevance 
is weak, in the former case. 

2.6 Oxygen abundance distribution in the SN halo 

The RaH07 subsample is considered to infer the EDOD 
related to the SN halo only for illustrative purposes, due 
to its incompleteness. The oxygen abundance range is 
-1.00 < [0/H] < 0.30 (LTE), -1.20 < [0/H] < 0.20 
(NLTE). 

The EDOD related to the RaH07 subsample is hsted 
in Tableiniand plotted in Fig.[3]both in presence (upper 
right) and in absence (lower right) of the LTE approx- 
imation. 

It can be seen that EDODs deduced from RaII07 and 
MaB08 subsamples show qualitative agreement. 

The EDOD related to the RaU07 subsample is also 
plotted in Fig. [3] for comparison, both in presence (up- 



per left) and in absence (lower left) of the LTE approx- 
imation. The related trend appears to be closer to its 
SN thick disk and SN thin disk counterparts, with re- 
spect to the SN halo. 

3 Inferred SN thick + thin disk oxygen 
abundance distribution 

Under the assumption of a universal initial mass func- 
tion (IMF) for a star generation, the SN thick + thin 
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Table 4 The empirical, diflerential oxygen abundance distribution (EDOD) in the solar neighbourhood (SN) thin disk, 
deduced from the RaN07 subsample (iV = 310) both in presence (LTE) and in absence (NLTE) of the local thermodynamic 
equilibrium approximation 
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Table 6 The empirical, differential oxygen abundance distribution (EDOD) in the solar neighbourhood (SN) halo, deduced 
from the RaH07 subsample {N = 28) both in presence (LTE) and in absence (NLTE) of the local thermodynamic equilibrium 
approximation 
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Fig. 2 The empirical, differential oxygen abundance distribution (EDOD) related to the (from top left in clockwise sense) 
MaKOS, MaN08, MaB08, subsample, and a combination of MaK08 and MaN08 subsamples after weighting by mass, for an 
assumed solar neighbourhood (SN) thick disk to thin disk mass ratio, Mkd/Mnd = 0.1 (diamonds) and 0.9 (triangles). The 
theoretical, differential oxygen abundance distribution (TDOD) related to inhomogeneous models of chemical evolution for 
the SN thick disk (case DB15), the SN thin disk (case DN81), the SN thick + thin disk (case KN4), and the bulge (case B) 
are plotted as full curves. The dotted curve on each panel reproduces (for comparison) the full curve on the other panel of 
the same row. The dashed curve plotted on the bottom right panel corresponds to a different bulge model (case Bl) 



disk EDOD may be expressed as (C07, COS): 



V' = log 

A=Fi/, = 

■ip^ = log 



Mkd A^Vkd , Mnd AiV] 



ND 
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log 
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A/nd AA^nd T (AA^nd)^/^ 
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KD 



Mr 



Ni 



ND 



-<5d) 



where Mkd and A/nd are the SN thick disk and the 
SN thin disk mass, Mr, = Mkd + A/nd, ^kd and 
TVnd the number of related sample objects, and ATVkd, 
AA'nd, the number of related sample objects within 
a selected metallicity bin, A[0/H], converted into A^ 
using Eq. ^ . For further details refer to the parent 
papers (C07, COS). 



In general, the disk is usually conceived as made of 
two main subsystems: the thick disk and the thin disk. 
Accordingly, the EDOD related to the SN disk depends, 
via Eqs. (O, on the SN thick to thin disk mass ratio, 
-^^kd/A/nd, which is poorly known at present. Values 
already quoted in literature s pan a wide r ange, from 
some perc ent (e.g.. Solmberg et al. ( 2007 )) to about 
unity (e.g., Fuhrmannl (|2008l )). or even indetermi nate in 



the se ns e that no di stinction can be made (e.g.. iNorris 



(5c) |l9S7t l: llvezic et all (jioOS) ). In addition, the global 



thick to thin disk mass ratio could be different from its 
SN counterpart, A/kdZ-^nd- 

The cases, A/kd/Mnd = 0.1, 0.3, 0.5, 0.7, 0.9, de- 
duced from RaK07 and RaN07 subsamples, are repre- 
sented in Fig. [Has triangles, diamonds, crosses, squares, 
asterisks, respectively, both in presence (LTE) and in 
absence (NLTE) of the LTE approximation. Error bars 
are plotted in left panels only, to allow comparison. 

The cases, A/kd/A/nd = 0.1 and 0.9, deduced from 
MaKOS and MaNOS subsamples, are represented in 
Fig. [2] (bottom left) as diamonds and triangles, respec- 
tively. 
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Fig. 3 The empirical, difFerential oxygen abundance distribution (EDOD) related to the RaU07 (left panels) and RaH07 
(right panels) subsample, both in presence (LTE, top panels) and in absence (NLTE, bottom panels) of the local ther- 
modynamic equilibrium approximation. The theoretical, differential oxygen abundance distribution (TDOD) related to 
inhomogeneous models of chemical evolution for the solar neighbourhood (SN) thick disk (case DB15), the SN thin disk 
(case DN), and the SN halo (case H) are plotted as full curves, dot-dashed curves, and long-dashed curves, respectively. 
Models are fitted to NLTE data. Bulge (case Bl) and SN halo (case HI) counterparts presented in an earlier attempt (C07) 
are also shown as dotted and dashed curves, respectively. 



An inspection of Figs. [2] (bottom left) and H shows 
that, in both cases, the uncertainty on the thick to thin 
disk mass ratio is comparable with the error box, except 
for the low-metallicity tail ([Fc/H]< —1) and/or poorly 
populated oxygen abundance bins, where a weak sta- 
tistical relevance is expected. 



4 Inhomogeneous, simple models 

In the light of inhomogeneous simple models of chem- 
ical evolution, a selected system is conceived as being 
structured into a number of discrete, entirely gaseous, 
identical regions, and a background of long-lived stars, 
stellar remnants, and gas inhibited from star forma- 
tion, which have been generated earlier. The evolu- 
tion occurs via a sequence of identical time steps. At 
the beginning of each step, star formation stochasti- 
cally takes place in a subclass of "active" regions, as 
described by simple homogeneous models, while the 
others remain "quiescent". At the end of each step, 
high-mass stars have died whereas low-mass stars have 



survived up until today, according to instantaneous re- 
cycling approximation. In addition, the enriched gas 
which remains from active regions is instantaneously 
mixed with the unenriched gas within quiescent regions, 
to form a new set of identical regions for the next step. 
For further details and complete formulation refer to 
earlier attemp t s (CO l, C07, COS) and parent papers 
( Mahnie et all (|l993h : COO). Due to reasons of simplic- 
ity, the current investigation shall be restricted to the 
special case of expected evolution, where the fraction 
of active regions is time independent. For the general 
case, refer to earlier attempts (COO, COl). 

Input parameters of inhomogeneous models consid- 
ered in the current attempt and parent papers (COO, 
COl, C07, COS) are: the number of steps during and af- 
ter the assembling phase; the IMF power-law exponent; 
the solar oxygen abundance; the minimum (beginning 
of the first step) and maximum (end of the last step) 
normalized oxygen abundance; the maximum normal- 
ized oxygen abundance at the end of the first step after 
the assembling phase; the true and the effective yield; 
the total number of regions at any step (if taken con- 



Oxygen abundance and chemical evolution 



11 



■s- 
<: 
z 



2; 
o 




LIE 



LIE 



-s- 
<i 

<i 

O 




NLTE 



NLTE 



-s- 
<i 
z 

<l 

bi 

o 




-s- 
<l 
z 

< 
o 



Fig. 4 The empirical, differential oxygen abundance distribution (EDOD) in the solar neighbourhood (SN) thick + thin 
disk, deduced from RaK07 and RaN07 subsamples, both in presence (LTE, top panels) and in absence (NLTE, bottom 
panels) of the local thermodynamic equilibrium approximation. Different symbols are related to different thick to thin disk 
mass ratios, Mkd/A/nd- Caption of symbols: Mkd/Mnd = 0.1 (triangles); 0.3 (diamonds); 0.5 (crosses); 0.7 (squares); 
0.9 (asterisks). Error bars are omitted in right panels to gain clarity. The theoretical, differential oxygen abundance 
distribution (TDOD) related to inhomogeneous models of chemical evolution for the SN thick + thin disk (cases KN), is 
plotted as dashed (case KN2) and full (case KN3) curves. Its counterpart, related to the SN thick disk (cases DB), is 
plotted as dotted (case DB105) curves. The dashed vertical bands on right panels corresp ond to [Fe/H1= — 1 a nd related 
uncertainties, deduced from Eq. @ with regard to a thick disk and a thin disk subsample (|Ramirez et al.l l|2007l )) 



stant during the evolution) or the mass of a region (if 
taken constant during the evolution); the total number 
of steps; as shown in Table[71 The more relevant output 
parameters are listed in Tables [H and [TOl A change 
in power-law IMF exponent affects only the mass frac- 
tion of a star generation, a, which remains locked up 
in long-lived stars and stellar remnants, and the lower 
stellar mass limit, m^{, while the remaining parameters 
are left unchanged. For further details see the parent 
papers (COl, COT, COS). 

Values of input parameters for different models re- 
lated to the thick disk (models DB), the thin disk (mod- 
els DN), the thick -I- thin disk (models KN), the halo 
(models H), and the bulge (models B) are listed in Ta- 
ble[7l Models Bl and HI come from earlier work (COT). 

Available samples, from which the EDOD has been 
deduced, are biased (towards low metallicities) and/or 
incomplete: for this reason, little meaning would be 
related to the best fit to the data. Different models 
considered here are mainly aimed in (i) analysing the 



effect of different input parameters on oxygen enrich- 
ment; (ii) providing acceptable fits to the data; (iii) 
investigating if a class of models with strictly universal 
IMF exist, which reproduce the EDOD determined for 
different populations. 

An equal number of steps has been chosen for both 
the assembling phase and the whole evolution, with re- 
gard to the halo, the bulge, and the thick disk , whic h 
formed withi n abo u t 1 Gyr (e.g., Prantzod ( 20081) ; 



Zoccali et ahl |2008l ); iMelendez et all (|2008l) ) On the 



other hand, a larger number of steps has been assumed 
for the thin disk, which is still undergoing chemical evo- 
lution. Accordingly, 4 episodes of star formation (2 
during the assembling phase) have been (somewhat ar- 
bitrarily) chosen for the halo, the bulge, and the thick 
disk, and 25 (5 during the assembling phase) for the 
thin disk. 

The initial and final normalized oxygen abundance 
have been chosen equal to the minimum and maximum 
value related to the most populated sample. With re- 
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Table 7 Values of input parameters related to the expected evolution of inhomogeneous simple models, for a number 
of different cases concerning the thick disk (models DB), the thin disk (models DN), the thick + thin disk (models KN), 
the halo (models H), and the bulge (models B). The expected evolution is independent of the total number of regions 
(COO). The assumed solar oxygen abundance is {Zo)q = 0.0056. The duration of evolution is needed only for specifying 
the age-metallicity relation (C07, COS). Different values of the fMF power-law exponent, p, leave the output parameters 
unchanged with the exception of the mass fraction of a star generation, a, which remains locked up in long-lived stars and 
stellar remnants, and the lower stellar mass limit, mmf, as shown in Tables[8l|9l and 1101 Symbol captions: Lc - number of 
steps during the assembling phase; La - number of steps during the whole evolution; ^min - minimum normalized oxygen 
abundance; 0max - maximum normalized oxygen abundance; ^^ax - maximum normalized oxygen abundance during the 
(Lc -|- l)-th step; a - slope of the empirical differential oxygen abundance distribution of the adjoint simple homogeneous 
model, see Eq. ([6l; p - true yield; p" - effective yield. The high number of digits in ^Jnax is necessary to ensure the same 
IMF during and after the assembling phase. For further details refer to the parent papers (COO, COl, C07) 



case 


Lc 


La 


0min 


0max 




-ap"/p 


P"/P 


DB12 


2 


4 


0.160 


2.5 


2.325653322 


0.48 


0.40 


DB13 


2 


4 


0.160 


2.5 


2.303503529 


0.52 


0.40 


DB14 


2 


4 


0.160 


2.5 


2.284610521 


0.56 


0.40 


DB15 


2 


4 


0.160 


2.5 


2.268320468 


0.60 


0.40 


DN 


5 


25 


0.400 


3.2 


2.01018844 


0.60 


0.40 


DN4 


5 


25 


0.400 


3.2 


2.38931821 


0.24 


0.60 


DN5 


5 


25 


0.400 


3.2 


2.13331046 


0.30 


0.60 


DN6 


5 


25 


0.400 


3.2 


1.96266176 


0.36 


0.60 


DN7 


5 


25 


0.400 


3.2 


1.84079119 


0.42 


0.60 


KNl 


5 


25 


0.160 


3.2 


2.09520539 


0.46 


0.40 


KN2 


5 


25 


0.160 


3.2 


1.83493684 


0.60 


0.40 


KN3 


5 


25 


0.160 


3.2 


1.59805348 


0.46 


0.60 


KN4 


5 


25 


0.160 


2.5 


1.90667713 


0.46 


0.40 


KN5 


5 


25 


0.160 


2.5 


1.420631883 


0.46 


0.60 


Bl 


2 


4 


0.200 


5.5 


4.44130682 


0.40 


0.70 


HI 


2 


4 


0.001 


1.0 


0.9834230800463403 


3.25 


0.04 


DB105 


2 


4 


0.160 


2.5 


2.367027589 


0.42 


0.40 


DN81 


5 


25 


0.200 


2.5 


1.51408616 


0.42 


0.60 


KN6 


5 


25 


0.160 


2.5 


2.01304781 


0.42 


0.40 


B 


2 


4 


0.200 


5.5 


4.43357173 


0.42 


0.70 


H 


2 


4 


0.100 


1.795 


1.7907654704 


0.42 


0.40 
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gard to the thin and the thick + thin disk, the above 
mentioned values are changed in some models (DN81, 
KN4, KN5, KN6) for improving the fit to the data. 
The maximum normalized oxygen abundance during 
the {Lc+l)-th step, (/)^ax: is fixed by the boundary 
condition of equal IMF during and after the assembling 
phase. The effective yield to the true yield ratio, p"/p, 
and the slope of the TDOD of the adjoint simple homo- 
geneous model, a, have been selected within a narrow 
range, outside which no acceptable fit to the data has 
been found. More specifically, the following considera- 
tions can be made. 

The lower and upper normalized oxygen abundance 
limit, (/'mill and 0max, have to be conceived as bound- 
ary conditions with respect to the subsystem under 
consideration. Assuming, say, 0niin = 0, would be 
meaningless in that (i) no star with zero metallicity 
has still been detected, and (ii) it would imply a G- 
dwarf problem in absence of ad hoc assumptions. Val- 
ues ((/>min, <?!'max) = (0.16,2.5) for the SN thick disk, 
(0.4,3.2) for the SN thin disk, (0.16,3.2) for the SN 
thick + thin disk, (0.001, 1.0) for the SN halo, (0.2, 5.5) 
for the bulge are taken from the most populated sam- 
ple from which the EDOD has been deduced. Values 
(0.2,2.5) for the SN thin disk and (0.16,2.5) for the 
SN thick 4- thin disk are obtained by removing a single 
oxygen overabundant star (HIP68184) from the RaN07 
subsample and, in the former case, assuming the low- 
metallicity tail of SN (still undetected) thin disk stars 
extends down to (j) = 0.2. Values (0.1,1.795) for the 
SN inner halo are taken from the incomplete RaHOT 
subsample, which makes the associated model poorly 
representative, even if conceptually interesting. 

The maximum normalized oxygen abundance dur- 
ing the {Lc + l)-th step, (j)'^;^^, defines the IMF during 
(first Lq steps) and after (subsequent La — Lc steps) 
the assembling phase. The high number of digits en- 
sures the same IMF during the whole evolution. It can 
be seen that the same number of steps during and af- 
ter the assembling phase, has been assigned to both 
the SN thick disk, the SN halo, and the bulge, on one 
hand, and to both the SN thin and the SN thick -I- thin 
disk, on the other hand. Strictly speaking, it would not 
imply a similar for mation timescale. as su ggested by re- 
cent findings (e.g., Melendez et al. ( 20081 )): to this aim, 
an equal duration would also be assigned to each step, 
which is relevant in dealing with the age-metallicity re- 
lation, but leaves the TDOD unchanged. 

The remaining input parameters are the effective to 
true yield ratio, p" /p, and the product: 



p" ^ (Zo)q 1 
p In 10 p 



(6) 



which is inversely proportional to the true yield. For 
further details refer to parent papers (COO, COl, C07). 

Values of output parameters for models shown in Ta- 
ble [7] are listed in Tables [8] and [H respectively. Ad- 
ditional cases where gas in neither inhibited from, nor 
enhanced in, star formation, are listed in Table [TUl The 
meaning of the parameters is explained in Appendix [XI 
For further details refer to parent papers (COO, COl, 
C07). 

The TDOD related to models DB15 and DN is com- 
pared in Fig.m to the corresponding EDOD with re- 
gard to SN thick (RaK07 subsample, left panels) and 
SN thin (RaN07 subsample, right panels) disk stars, 
respectively, both in presence (top panels) and in ab- 
sence (bottom panels) of the LTE approximation. On 
each panel, the fit to the data is plotted as a full curve, 
and its counterpart related to the right or left panel 
as a dashed line for comparison. Models are fitted to 
NLTE data. 

The TDOD related to models DB15, DN81, KN4, 
and B, is compared in Fig. [2] to the corresponding 
EDOD with regard to SN thick disk (MaK08 subsam- 
ple, top left), SN thin disk (MaN08 subsample, top 
right), SN thick + thin disk [MaK08 -I- MaN08 subsam- 
ple with assumed A/kdZ-^nd = 0.1 (diamonds) and 0.9 
(triangles), bottom left], and bulge (MaB08 subsample, 
bottom right) stars, respectively. On each panel, the fit 
to the data is plotted as a full curve, and its counter- 
part related to the right or left panel as a dotted line 
for comparison. Model Bl is represented by the dashed 
curve (bottom right). 

The TDOD related to models DB15 (fuU), DN (dot- 
dashed), Bl (dotted), and HI (dashed), is compared in 
Fig. [3] to the corresponding EDOD with regard to SN 
uncertain population (RaU07 subsample, left panels) 
and SN halo (RaII07 subsample, right panels) stars, re- 
spectively, both in presence (top panels) and in absence 
(bottom panels) of the LTE approximation. 

The TDOD related to models KN2 (dashed) and 
KN3 (full) is compared in Fig. [3] to the corresponding 
EDOD, with regard to SN thick -I- thin disk (RaK07 -I- 
RaN07 subsample) stars for assumed SN thick to thin 
disk mass ratio, A/kd/A/nd = 0.1 (triangles), 0.3 (dia- 
monds), 0.5 (crosses), 0.7 (squares), and 0.9 (asterisks). 
Both error bars and TDODs are omitted in right pan- 
els to gain clarity. Model DB105, related to the SN 
thick disk, is also represented as a dotted line. The 
dashed vertical bands correspond to [Fe/II]= —1 and 
related uncertainties, according to Eq. with regard 
to a thick disk (left) and a thin disk (right) subsam- 
ple (|Ramirez et all |2007t )). Models are fitted to NLTE 
data. 

The TDOD related to models DB12, DB13, DB14, 
DB15, DB105 (from up to down along the vertical line, 
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Table 8 Values of output parameters related to the expected evolution of inhomogeneous simple models, for four different 
cases concerning the thick disk (models DB), and five concerning the thin disk (models DN), respectively. The indices, 
2.9 and 2.35, denote values related to the power-law IMF exponent, p, in computing the corresponding quantities. For the 
parameter, tpi, upper and lower values are calculated as in an earlier attempt (C07) by use of Eqs. (40) and (41) therein, 
respectively. The effective yield, p' , is related to inhomogenities in oxygen abundance due to the presence of active and 
quiescent regions, whereas oxygen is uniformly distributed within active regions. The lower part of the Table (last 6 rows) 
is related to models with inhibited star formation. Parameters not reported therein have the same value as in the upper 
part, with the exception of p, a, and m^f = mmf/m©, which are listed in Table [TOl together with other parameters not 
appearing here. The effective yield, p" , due to the presence of (inhibiting star formation) gas within active regions, is listed 
as p in the upper part of the table. The effective yield, p', due to the presence of both (inhibiting star formation) gas within 
active regions, and (precluding star formation) gas within quiescent regions, is listed with the same notation in the upper 
part of the table. The index, R, denotes a generic active region. The mean oxygen abundance (normalized to the solar 
value) of stars at the end of evolution is denoted as (j) 

I LTE I NLTE 



DB12 DB13 DB14 DB15 DN DN4 DNS DN6 



Mr 


1 


.3243 


E- 


1 


1, 


.2295 


E- 


1 


1, 


.1258 


E- 


1 


1, 


.0304 


E- 


1 


1, 


.6663 


E- 


1 


3, 


.7477 


E- 


1 


3 


.6662 


E- 


-1 


3, 


.5864 E- 


q 


6 


.1771 


E- 


1 


5, 


.5533 


E- 


1 


4, 


.9941 


E- 


1 


4, 


.4924 


E- 


1 


8, 


.0550 


E- 


1 


9, 


.6146 


E- 


1 


9 


.3741 


E- 


1 


9, 


.1396 E- 


X 


4, 


.41S6 


E- 


1 


5, 


.0700 


E- 


1 


5, 


.6410 


E- 


1 


6, 


.1403 


E- 


1 


2, 


.3339 


E- 


1 


6, 


.1635 


E- 


-2 


9 


.8821 


E-^ 


2 


1, 


.3416 E- 




9 


.0478 


E- 


1 


8, 


.3518 


E- 


1 


7, 


.7553 


E- 


1 


7, 


.2382 


E- 


1 


7 


.2382 


E- 


1 


1, 


.8096 


E- 


-0 


1 


.4476 


E- 





1, 


.2064 E- 




3 


.6191 


E- 


1 


3, 


.3407 


E- 


1 


3, 


.1021 


E- 


1 


2, 


.8953 


E- 


1 


2. 


.8953 


E- 


1 


1, 


.0857 


E- 


-0 


8, 


.6859 


E- 


1 


7, 


.2382 E- 




1 


.4SS9 


E- 


1 


9, 


.5108 


E- 


2 


6, 


.2204 


E- 


2 


4, 


.0731 


E- 


-2 


4, 


.4842 


E- 


3 


3, 


.7439 


E- 


1 


1, 


.9872 


E- 


1 


1, 


.0548 E- 


0:2.9 


6 


.9S0S 


E- 


1 


7, 


.1175 


E- 


1 


7, 


.2671 


E- 


1 


7, 


.4020 


E- 


1 


7, 


.4020 


E- 


1 


5, 


.3263 


E- 


1 


5, 


.8755 


E- 


1 


6, 


.3092 E- 


^2.35 


8, 


.6951 


E- 


1 


8, 


.7833 


E- 


1 


8, 


.8803 


E- 


1 


8, 


.9281 


E- 


1 


8, 


.9281 


E- 


1 


7, 


.6914 


E- 


1 


8, 


.0638 


E- 


-1 


8, 


.3327 E- 


(Wmf)2.9 


4, 


.0227 


E- 


1 


3, 


.7813 


E- 


1 


3, 


.5660 


E- 


1 


3 


.3728 


E- 


1 


3 


.3728 


E- 


1 


6, 


.4228 


E- 


1 


5, 


.6022 


E- 


1 


4, 


.9602 E- 


(minf)2.35 


1 


.130S 


E- 


2 


9, 


.3463 


E- 


3 


7, 


.8194 


E- 


3 


6 


.6108 


E- 


3 


6 


.6108 


E- 


3 


5, 


.1835 


E- 


-2 


3 


.2618 


E- 


-2 


2, 


.1898 E- 




1 


.743S 


E- 


1 


1, 


.9650 


E- 


1 


2, 


.1539 


E- 


1 


2, 


.3168 


E- 


1 


6 


.2622 


E- 


2 


4, 


.2667 


E- 


2 


5, 


.6142 


E- 


2 


6, 


.5123 E- 




1 


.8170 


E- 





1, 


.7505 


E- 





1, 


.6938 


E- 





1, 


.6450 


E- 





1, 


.2971 


E- 





1, 


.7760 


E- 





1 


.4526 


E- 





1, 


.2370 E- 


-^0 


2, 


.4322 


E- 


1 


1, 


.7337 


E- 


1 


1, 


.1035 


E- 


1 


5, 


.3384 


E- 


-2 


4, 


.8961 


E- 


-2 


1, 


.2641 


E- 





1 


.0696 


E- 





9, 


.0546 E- 


-^-1(40) 


2, 


.8439 


E- 


1 


2, 


.2346 


E- 


1 


1, 


.6926 


E- 


1 


1, 


.2104 


E- 


1 


5, 


.0826 


E- 


1 


1, 


.2692 


E- 


-0 


1 


.0780 


E- 





9, 


.1713 E- 


-V'i(41) 


2, 


.8S06 


E- 


1 


2, 


.2446 


E- 


1 


1, 


.7066 


E- 


1 


1, 


.2289 


E- 


1 


5, 


.0839 


E- 


1 


1, 


.2692 


E- 


-0 


1 


.0780 


E- 





9, 


.1718 E- 




1 


.4286 


E- 


1 


2, 


.3810 


E- 


1 


3, 


.3333 


E- 


1 


4, 


.2857 


E- 


1 


4, 


.2857 


E- 


1 


-4, 


.2857 


E- 


T 


-2 


.8571 


E- 


1 


-1, 


.4286 E- 




7 


.5755 


E- 


1 


7, 


.0839 


E- 


1 


6, 


.6557 


E- 


1 


6 


.2787 


E- 


1 


5, 


.8336 


E- 


1 


1, 


.0942 


E- 





8, 


.8673 


E- 


1 


7, 


.4825 E- 




1 


.0822 


E- 


1 


1, 


.6866 


E- 


1 


2, 


.2186 


E- 


1 


2, 


.6909 


E- 


1 


2, 


.5001 


E- 


1 


-4, 


.6892 


E- 


1 


-2 


.5335 


E- 


1 


-1, 


.0689 E- 




7 


.8307 


E- 


1 


7 


.4978 


E- 


1 


7, 


.2065 


E- 


1 


6 


.9485 


E- 


1 


8, 


.6448 


E- 


1 


1, 


.1450 


E- 


-0 


1, 


.0068 


E- 





9, 


.1462 E- 




7. 


.4761 


E- 


1 


7, 


.3087 


E- 


1 


7, 


.0335 


E- 


1 


6 


.7149 


E- 


1 


6, 


.9686 


E- 


1 


1, 


.0948 


E- 


-0 


1 


.1218 


E- 





1, 


.0436 E- 


Df 


1 


.0680 


E- 


1 


1, 


,7402 


E- 


1 


2, 


.3445 


E- 


1 


2, 


.8778 


E- 


1 


2, 


.9865 


E- 


1 


-4, 


.6921 


E- 


1 


-3 


.2051 


E- 


1 


-1, 


.4909 E- 
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Table 9 Values of output parameters related to the expected evolution of inhomogeneous simple models, for five different 
cases concerning the thick + thin disk (models KN). Additional cases related to the bulge (model Bf) and the halo (model 
Iff), taken from an earlier attempt (C07), but related to a different reference model with respect to inhibition from and 
enhancement in star formation, are listed for comparison. Captions as in Table |8] 



KNl KN2 KN3 KN4 KN5 Bl HI 





1 


7520 


E- 


1 


1 


6240 


E- 


1 


3 


4073 


E- 


1 


1 


8550 


E- 


1 


3 


5544 


E- 


1 


1 


4139 


E- 


1 


8 


2201 


E- 


4 


5 


8 


5730 


E- 


1 


7 


8025 


E- 


1 


8 


6171 


E- 


1 


9 


2064 


E- 


1 


9 


0458 


E- 


1 


2 


4833 


E- 


1 


4 


4991 


E- 


2 


X 


1 


7302 


E- 


1 


2 


6236 


E- 


1 


2 


0977 


E- 


1 


9 


7438 


E- 


2 


1 


4804 


E- 


1 


8 


7545 


E- 


1 


9 


5579 


E- 


1 


vI{Zo)q 


9 


4412 


E- 


1 


7 


2382 


E- 


1 


9 


4412 


E- 


1 


9 


4412 


E- 


1 


9 


4412 


E- 


1 


1 


0857 


E- 





1 


3363 


E- 


1 


p'I{Zo)q 


3 


7765 


E- 


1 


2 


8953 


E~ 


1 


5 


6647 


E- 


1 


3 


7765 


E- 


1 


5 


6647 


E- 


1 


7 


6001 


E- 


1 


5 


3452 


E- 


3 




2 


1295 


E- 


2 


2 


0221 


E- 


3 


2 


4210 


E- 


2 


1 


2653 


E- 


1 


8 


1501 


E- 


2 


3 


8030 


E- 


3 


4 


0973 


E- 


6 


02.9 


6 


8596 


E- 


1 


7 


4020 


E- 


1 


6 


8596 


E- 


1 


6 


8596 


E- 


1 


6 


8596 


E- 


1 


6 


5510 


E- 


1 


9 


3915 


E- 


1 


Q;2.35 


8 


6460 


E- 


1 


8 


9281 


E- 


1 


8 


6460 


E- 


1 


8 


6460 


E- 


1 


8 


6460 


E- 


1 


8 


4739 


E- 


1 


9 


7832 


E- 


1 


(TOmf)2.9 


4 


1547 


E- 


1 


3 


3728 


E- 


1 


4 


1547 


E- 


1 


4 


1547 


E- 


1 


4 


1547 


E- 


1 


4 


6050 


E- 


1 


6 


7781 


E- 


2 


(TOmf)2.35 


1 


2495 


E- 


2 


6 


6108 


E- 


3 


1 


2495 


E- 


2 


1 


2495 


E- 


2 


1 


2495 


E- 


2 


1 


7263 


E- 


2 


7 


5651 


E- 


5 




5 


8147 


E- 


2 


7 


1945 


E- 


2 


8 


4313 


E- 


2 


3 


1228 


E- 


2 


5 


6809 


E- 


2 


1 


0587 


E- 





1 


6577 


E- 


2 




1 


6445 


E- 





1 


3157 


E- 





1 


0165 


E- 





1 


5905 


E- 





9 


7659 


E- 


1 


2 


1239 


E- 





9 


4927 


E- 


1 


-■00 


7 


2759 


E- 


1 


4 


3986 


E- 


1 


6 


4265 


E- 


1 


9 


2754 


E- 


1 


7 


6772 


E- 


1 


9 


1837 


E- 


2 


8 


5447 


E- 


1 


-V'i(40) 


7 


4089 


E- 


1 


4 


6123 


E- 


1 


6 


6189 


E- 


1 


9 


3471 


E- 


1 


7 


8072 


E- 


1 


2 


8650 


E- 


1 


8 


2781 


E- 


1 


-V^i(41) 


7 


4096 


E- 


1 


4 


6141 


E- 


1 


6 


6204 


E- 


1 


9 


3473 


E- 


1 


7 


8079 


E- 


1 


3 


0358 


E- 


1 


8 


2753 


E- 


1 


K 


9 


5238 


E- 


2 


4 


2857 


E- 


1 


9 


5238 


E- 


2 


9 


5238 


E- 


2 


9 


5238 


E- 


2 


-4 


7619 


E- 


2 


6 


7381 


E- 







7 


5307 


E- 


1 


5 


8632 


E- 


1 


6 


0194 


E- 


1 


7 


4367 


E- 


1 


5 


8851 


E- 


1 


9 


0154 


E- 


1 


1 


2912 


E- 


1 


D'k 


7 


1721 


E- 


2 


2 


5128 


E- 


1 


5 


7327 


E- 


2 


7 


0826 


E- 


2 


5 


6048 


E- 


2 


-4 


2930 


E- 


2 


8 


7005 


E- 


1 


<!> 


7 


5480 


E- 


1 


6 


2873 


E- 


1 


5 


7876 


E- 


1 


7 


4187 


E- 


1 


5 


6558 


E- 


1 


9 


3597 


E- 


1 


1 


3385 


E- 


1 




8 


9360 


E- 


1 


6 


9858 


E- 


1 


8 


9094 


E- 


1 


7 


9751 


E- 


1 


8 


3863 


E- 


1 


1 


0460 


E- 





1 


2923 


E- 


1 


Df 


8 


5105 


E- 


2 


2 


9939 


E- 


1 


8 


4851 


E- 


2 


7 


5954 


E- 


2 


7 


9869 


E- 


2 


-4 


9810 


E- 


2 


8 


7077 


E- 


1 
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Table 10 Values of output parameters related to the expected evolution of inhomogeneous simple models, for five different 
cases concerning the thick disk (model DB105), the thin disk (model DN81), the thick + thin disk (model KN6), the bulge 
(model B), and the inner halo (model H). The stellar initial mass function (IMF) is strictly universal, in the sense that 
neither inhibition from, nor enhancement in, star formation, appear in one model with respect to the other. Captions as in 
Table E] 

DB105 DN81 KN6 B H 
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-0 


p'I{Zo)q 


4, 


.1361 


E- 


1 


6 


.2042 


E- 


1 


4, 


.1361 


E- 


1 


7, 


.2382 


E- 


1 


4, 
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E- 
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E- 
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8, 
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E- 
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8, 
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8, 
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8, 
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9, 
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.0000 


E- 





0, 
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Fig. 5 The empirical, differential oxygen abundance distribution (EDOD) related to the RaKOT (top left), RaN07 (top 
right and bottom left) and RaK07 + RaN07 (bottom right) subsample, for assumed SN thick to thin disk mass ratio, 
Mkd/Mnd = 0.1 (diamonds) and 0.9 (triangles), in absence of the local thermodynamic equilibrium (LTE) approximation. 
The theoretical, differential oxygen abundance distribution (TDOD) related to inhomogeneous models of chemical evolution, 
is plotted for the following models; DB12, DB13, DB14, DB15, DB105 (from up to down along the vertical line, <j) = 2.0, top 
left); DN4, DN5, DN6, DN7, DN (from up to down along the vertical line, = 2.5, top right); DN, DN81 (from up to down 
along the vertical line, 4> = 2.4, bottom left); KN2, KN4, KN6 (from down to up along the vertical line, 4> ~ 2.0, bottom 
right). The dashed vertical bands correspond to [Fe/H]= —1 and related uncertainties, deduced fro m Eg. (l4ll with regard t o 
a thick disk (upper left and lower right) and a thin disk (upper right and bottom panels) subsample (jRamirez et al.l (|2007l )) 



(j) = 2.0, top left); DN4, DN5, DN6, DN7, DN (from 
up to down along the vertical line, (j) — 2.5, top right); 
DN, DN81 (from up to down along the vertical line, 
<p = 2.4, bottom left); KN2, KN4, KN6 (from down to 
up along the vertical line, (j) — 2.0, bottom right); is 
compared in Fig. [5] to the corresponding EDOD with 
regard to SN thick disk (RaK07 subsample, top left), 
SN thin disk (RaN07 subsample, top right and bot- 
tom left), and SN thick + thin disk [RaKOT + RaNOT 
subsample with assumed A/kdZ-^nd = 0.1 (diamonds) 
and 0.9 (triangles), bottom right] stars, respectively, in 
absence of the LTE approximation. 



5 Discussion 

5.1 Fitting empirical differential oxygen abundance 
distributions 

Acceptable fits to the EDOD related to SN thin or thick 
+ thin disk stars, must satisfy the boundary condition 



of a nonnegligible prese nt-day gas mass fracti o n, [if = 
Mg/M = 0.1-0.3 (e.g., IPrantzos and AubertI (|l995l) ). 
Accordingly, related models (DN and KN, respectively) 
must be disregarded if (i) the EDOD is poorly fitted 
and/or (ii) /i/ lies outside the above mentioned range. 

With regard to the SN thin disk, an inspection of 
Figs.IUISliniand Table [5] shows that the sole acceptable 
model (DN81) needs lowered initial and final oxygen 
abundance, close to thick disk values. More specifi- 
cally, fitting the data implies the existence of SN thin 
disk stars with normalized oxygen abundance down to 
(j) « 0.2 and, in addition, the removal of a single oxy- 
gen overabundant ([0/H]nlte = 0.48 or <^ w 3) star 
(HIP68184) from the sample, which could be due to 
peculiar evolutionary effects and/or systematic errors. 
In the light of the model, the following conclusion holds: 
if the thick and the thin disk evolved as separate en- 
tities, the chemical evolution spanned along a similar 
oxygen abundance range. 
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With regard to the SN thick + thin disk, an inspec- 
tion of Fig. [5] and Table [9] shows that the two accept- 
able models (KN4 and KN6) need lowered final oxygen 
abundance, close to thick disk values. More specifically, 
fitting the data implies the removal of the above men- 
tioned, oxygen overabundant star, from the sample. In 
the light of the model, the following conclusion holds: 
if the thick and the thin disk evolved as a single entity, 
the chemical evolution spanned along a similar oxygen 
abundance range with respect to the thick disk consid- 
ered alone. Models where the thick and the thin disk 
evolve as a whole, are consistent with the existence of 
an evolutionary fi nk between the two disk components 
(Haywood ( 2008h ). It is worth remembering that the 
above quoted paper is referring to a kinematic evolution 
or migration between thin and thick disks, as opposed 
to some form of Galactic chemical evolution. 

With regard to the bulge, an inspection of Figs. [21 [3] 
and Tables [9l [TOl shows that both the models considered 
(B and Bl) provide acceptable fits to bulge and SN halo 
data, with negligible amount of gas left at the end of 
evolution, as expected. In the light of the model, the 
following conclusion holds: oxygen production within 
the bulge and the SN halo underwent a similar history 
for the common abundance range. This result is highly 
speculative, due to the incompleteness o f the halo sam 



pie fro m which the data were extracted ([Ramirez et al 
( 2007l )'l. 



With regard to the SN halo, an inspection of Fig.[3| 
and Tables [^ [TUl shows that a fraction of low-metallicity 
data is fitted by a single model (HI), while the remain- 
ing part is fitted by a second model (H). Accordingly, 
the (incom plete) halo sample fro m which the data were 
extracted ([Ramirez et al.[ (j2007[ )) has to be conceived 
as made of two distinct classes of objects, namely (i) 
stars from a gas reservoir from which both the inner 
halo, the bulge, and the thick disk originated, and (ii) 
stars from a gas reservoir from which the outer halo 
originated. In the light of the model, the following con- 
clusion holds: oxygen production within the SN inner 
halo, the bulge, and the thick disk underwent a sim- 
ilar history for the common abundance ranges, while 
the chemical evolution of the SN outer halo was differ- 
ent. In fact, input parameters related to the outer halo 
model (HI) are markedly different from the remaining 
ones, as shown in Table[7l A stellar dichotomy between 
the inner and the out er halo has recently been detected 
( Carollo et aP ^m^). 



The assumption of universal IMF implies gas inhibi- 
tion from, or enhancement in, star formation, with re- 
spect to a reference thin disk model (C07, COS), which 
has been chosen to be DN81. Then a legitimate ques- 
tion is about the existence of models which fit to an 



acceptable extent the EDOD related to the SN thick 
disk, SN thin disk, SN inner halo, bulge, and where, 
in addition, gas is neither inhibited from, nor enhanced 
in, star formation. The answer is positive, as shown in 
Table [TOl and Figs.[ll [3 and \B 

For a specified universal IMF, the key parameter ap- 
pears to be the probability of a region being active, 
X- The chemical evolution of the SN thick, thin, and 
thick + thin disk, appear to be characterized by low 
values, X = 0.32, 0.13, 0.07, respectively, while a high 
value, X = 0.89, is related to the chemical evolution 
of the bulge, as shown in Table [TUl In addition, the 
chemical evolution of the SN inner halo is character- 
ized by X = 0.013, a value close to one tenth its SN 
thin disk counterpart. The corresponding true oxygen 
yield, p = 0.0058, is consistent wi th jp = 0.0 085 T 0.0035 
deduced from observations (e.g., Prantzoj (fl994) ). 

In conclusion, a single inhomogeneous model of 
chemical evolution with universal IMF and gas neither 
inhibited from, nor enhanced in, star formation, pro- 
vides acceptable fits to the EDOD related to both the 
SN inner halo, the SN thick disk, the SN thin disk, the 
SN thick -I- thin disk, and the bulge. In particular, the 
main part of objects of uncertain population could be- 
long to the thin disk, as shown in Fig. [31 Conversely, an 
inhomogeneous model with considerable amount of gas 
inhibited from star formation (HI), used in an earlier 
attempt (C07), provides an acceptable fit to the outer 
halo. 

With regard to the SN thick disk, the SN thin disk, 
and the bulge, both poor and rich samples are fitted 
by the same models. With regard to the SN halo, 
the existence of two different classes of objects in the 
(incomplete) related sample, has to be advocated: a 
low-metallicity tail belonging to the outer halo, and 
a higher-metallicity body belonging to the inner halo, 
each fitted by a different model, as shown in Fig.[3l A 
possible explanation shall be given below. 

5.2 A physical interpretation of the model 

The EDODs related to MaK08, MaN08, MaK08 + 
MaN08, MaB08, RaK07, RaN07, RaK07 + RaN07 sub- 
samples, are fitted to an acceptable extent by models 
DB105, DN81, KN6, B, fisted in Table [TOl In addi- 
tion, the EDOD related to the RaH07 subsample is fit- 
ted to an acceptable extent by model H listed in Table 
[TOl provided the low-metallicity tail is not considered. 
Present-day gas mass fraction and oxygen abundance 
related to the SN disk may be expressed as weighted 
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means: 

M 



Kg 



iZo)f 



Mk 

Mk + A'/n 
(Mo)g , 



Mn 

AfN^Ng 



A^k + A/n 

_ (Afo)Kg + (Afo)Ng 

ATg AfKg + AfNg 

AfKg(^o)Kg A//Ng(^o)Ng 



Ml 



Ng 



Mv 



Ml 



Ng 



(7) 



(8) 



where Afg and Af are the present gas and total mass, 
and the indices, K and N, denote the thick and the thin 
disk, respectively. For models DB105 and DN81, the 
reference values, Mk/M-^ = 0.1, 1, yield /i/ = 0.13746, 
0.19998, respectively, to be compared to /i/ = 0.21244 
related to model KN6, where the thick and thin disk 
evolve as a single system. 

Models listed in Table (TU] are related to a strictly 
universal IMF, in the sense that both the power-law 
exponent and the upper and lower stellar mass limit 
coincide in all cases and, in addition, gas is neither 
inhibited from, nor enhanced in, star formation, con- 
cerning one model with respect to the other. This is 
a stronger condition in comparison with earlier results 
(COl, COT, COS) based on SN disk samples biased to- 
wards high-mctallicity stars. 

For models characterized by a strictly universal IMF, 
the probability of a region being active, x, is highly 
relevant. An inspection of Table (TU] shows that x — 
0.01268, 0.07405, for the SN inner halo and the SN thick 
+ thin disk, respectively, increasing up to x = 0.12549, 
0.32460, 0.88717, for the SN thin disk, the SN thick 
disk, and the bulge, respectively. A physical interpre- 
tation of the above results, in the light of the model, 
is that the probability of a region being active mainly 
depends on the environment. 

More specifically, the star formation rate within an 
active region is independent of the environment in the 
case under discussion. On the contrary, a different 
global efficiency (which is proportional to the proba- 
bility of a region being active) takes place in different 
environments. Using the values of x listed in Table 
[TOl the following global efficiency ratios are calculated: 
Xb/xh « 70; xb/xkn ~ 12; xb/xn « 7; xb/xk « 3; 
where B, H, K, N, KN denote the bulge, the SN halo, 
the SN thick disk, the SN thin disk, the SN thick + 
thin disk, respectively. 

Star formation is triggered by tidal interaction 
(including merging, which can be conceived as the 
strongest tidal interaction). Accordingly, lower x val- 
ues are expected in environments where tidal effects are 
poorly experienced, such as the inner halo, the thick 
+ thin disk, and the thin disk. On the other hand. 



larger x values are expected in environments where 
considerable tidal effects take place, such as the thick 
disk and the bu l ge. A ccording to recent findings (e.g., 
Melendez et al.l (|2008[ )') the thick disk and the bulge 
formation timescales were similar, of the order of sev- 
eral tenths of Gyr. The density is higher at the end of 
contraction, which could maximize tidal effects between 
different regions, and then produce a larger number of 
active regions with respect to environments where the 
density is considerably lower (inner halo) or the chemi- 
cal evolution timescale considerably longer (thin disk) . 

5.3 Universal IMF and G-dwarf problem 

A different IMF in bulge and disk components has 
been advocated for explaining the metallicity distri- 
bution in both the M ilky Way and Andromeda bulge 
iBallero et al.l (|2007bh ). contrary to what has been 
deduced from recent observations performed on SN 
thick disk, SN thin di sk, and Milky Way bulge stars 
( Melendez et al.l ( 2008 )). As pointed out in an earlier 
attempt (C08), a large fraction of Fe production de- 
pends on type la supernovae which, in turn, are thought 
to be born from an accreting white dwarf progenitor in 
a binary system, including a red giant outside its Roche 
lobe. Accordingly, different IMFs in different environ- 
ments could be simulated by different amounts of bi- 
nary star fractions (in particular, Snia progenitors). In 
this view, denser environments where objects with op- 
posite velocity vectors are close one to the other (e.g., 
the assembling bulge, the assembling thick disk) should 
be expected to exhibit a smaller binary star fraction (in 
particular, SnIa progenitors), with respect to less dense 
environments where objects with parallel velocity vec- 
tors are close one to the other (e.g., the assembled thin 
disk). 

Regions undergoing an early infall phase (of pri- 
mordial composition), with the smaller regions hav- 
ing a more important infall, have been advocated 
for explaining the field star metallicity distribution in 
both the halo and dwarf spheroidal satellites of the 
Milky Wa y , where a G-dwarf problem seems to occur 
|Prantzosl (|1994) 1. Though a G-dwarf problem (i.e. 
overabundance of low-metallicity, long-lived stars pre- 
dicted by the Simple model with respect to observa- 
tions) has been detected both in bulge-domina t ed an d 



lenrv and Wort 



lev 



1999h l 



disk-dominated galaxies 

and is probably universal ( Worthev et al.1 ( 19961 )). still 
it ma y be explained by a l arge n umber of different ways 
(e.g.. |Pagel and Patchettl()l975h ) including, among oth- 
ers, unprocessed gas inflow. On the other hand, a 
prompt initial enrichment from pop. Ill stars cannot 
be overlooked and, in particular, the idea of a coeval 
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formation of p op. Ill and pop. II stars cannot be dis- 
regarded (e.g., ISmith et al.l ([20081)). In this view, the 
G-dwarf problem in old stellar systems (e.g., the halo) 
could find a natural explanation, if still occurring within 
the framework of inhomogeneous models of chemical 
evolution. 

5.4 Stellar migration across the disk 

Though sample objects considered in the current pa- 
per, with the exception of bulge stars, are made of SN 
stars, still related informations provide valuable clues 
for understanding the formation process of the Galaxy. 
The assumption that the EDOD of the local disk is 
representative of the global disk, even if in contrast 
with an inside-out disk formation, can be considered 
as a useful zero-th order approximation. On the other 
hand, nearby stars older than about 0.2 Gyr come from 
birth sites whi ch span a large range in Galactocentric 



distances 



e-g-, 



Rocha-Pinto et al 



The orbital 



diffusion coefficient deduced from the observed increase 
of velocity dispersion with age implies that presently lo- 
cal stars have suffered a rms azimuthal drift from about 
2 kpc (for an age of 0.2 Gyr) to many Galactic orbits 
(for a n age of 1 Gyr ); for further details refer to earlier 
work dWielenl (|l977h 'l. Considerable but s maller, drift 



should occur also on the radial direction CWielen et al 
( 199(ih : lHavwoodl (|2008h V 

In this sense, the star formation rate inferred for 
nearby stars is a measure of the global Milky Way star 
form ation rate, at least at the sun Galactocentric ra- 
dius ( Rocha-Pinto et al. ( 2000f) ). acco rding to the est! 



mates of the diffusion coefficient (e.g.. iMeusinger et al 
(iiiil)). Star migration along the equatorial plane is 



mainly due to churning and blurring via interactions 
betwe en stars and spiral arms (Se llwood and Binney, 
2002; ISchonrich and BinnevI (|2009f l). Diffusion calcu- 
lations based on the increase of local velocity disper- 
sion with age, could be inaccurate if significant migra- 
tion occurs, as due to resonant scattering with tran- 
sient spiral arms, which makes local samples consid- 
erably affected by radial migration, according to A^- 
body -|- smooth pa rticle hydrodynamics simulations 



(Roskar et al 



(|2008h l. 

As outlined in an earlier attempt ( Wang and Silkl 
( 1993^ ). radial flows of the gas (possibly due to the 
observed spiral density waves) or different disk star for- 
mation histories (with time scales comparable to that 
of chemical evolution) between the inner and the outer 
parts of the Galaxy, may change the local abundances, 
but the overall abundance in the disk should not be 
affected as long as there is no gain o r loss of mate- 



rial in the disk. It can also be noticed (jWang and 



(il993;.) ) that the average oxygen abundance in the disk 
is roughly the solar value, and can be plausibly ex- 
plained by the standard Scalo IMF with lower and up- 
per star mass limit, fhmf = 0.1 and rfiMf = 60, re- 
spectively, or by a power-law IMF with exponent, p, 
within the range, 2.35 < p < 2.9, and the same fhmf 
and friMi- For further details refer to the parent paper 
( Wang and sii5 d 19931) 1. 

Inhomogeneous models can describe to some extent 
stellar migration across the disk, provided inflow (i.e. 
enhanchement of star formation) is considered. As out- 
lined in an earlier attempt (C07), inhomogeneous mod- 
els where star formation is enhanced imply star in- 
flow with same oxygen abundance as in the existing 
gas, which mimics (a special case of) stellar migration. 
Mass conservation no longer holds, and the effective 
gas mass fraction (normalized to the final mass), /Xoff, 
is expressed in terms of the actual gas mass fraction 
(normalized to the initial mass), /i, as: 



fJ-cS 



1 



(9) 



where Smig is the infiowed star mass fraction (normal- 
ized to the initial mass), to be interpreted as the net 
stellar migration within the volume under considera- 
tion, i.e. the difference between the fractional mass 
in stars which have been born elsewhere and migrated 
inside, and born inside and migrated elsewhere, respec- 
tively. 

In the framework of inhomogeneous models with en- 
hanced star formation, the fractional mass in migrated 
stars is proportional to the fractional mass in stars gen- 
erated in situ, Smig — (s = C(l ^ /^), and Eq. ^ may 
be rewritten as: 



M = 



(1 + pMcff 
1 + C^eff 



(10) 



where the effective gas mass fraction, fics, can be de- 
duced from the observations. For the solar neigh- 



bourhood, 0.1 ~ fj,cs ~ 0.3 (iPrantzos and Aubert 



( 19951 )). ( Ki I from high-resolution A^-body -I- smooth 
particle hydrodynamics simulations of disk formation 
iRoskar et al.l (|2008h ). and Eq. ^ yields 0.18 - ^ - 
0.46. 

On the other hand, model DN4 listed in Table [5] 
produces /i ~ 0.37, s = 1 — /i « 0.63 and and 
—Df — Smig ~ 0.47, which implies about three quarters 
of stars were born in situ, and the remaining quarter is 
the net effect of stellar migration. Accordingly, inhomo- 
geneous models with enhanced star formation provide 
an acceptable description of stellar migration even if, in 
the case considered, model DN4 fails in fitting the SN 
thin disk EDOD, as shown in Fig. [5] (top right panel). 
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5.5 Implications for the formation of the Galaxy 

According to a number of recent investigations, the for- 
mation of the Galaxy appears to be characterized by the 
following main features. 

(i) Thick disk stars had a chemical enrichment sim- 
ilar an d coeval to met a l-rich halo stars, [Fe/Hjsi 
-1.3 dProchaska et al.l tood^ ). -1.5 <[Fe/H]< 
2mih ). -1.8 <[Fe/H]< -0.5 



0.1 (Ramirez et 



al 



( Melendez et al.l (|2008[ )). 
(ii) Thick disk stars had a chemical enrichment sim- 
ilar a nd coeval to metal-poo r bulge stars, [Fe/II]~ 
- 0.4 dProchaska et al.l (l2000h). -1.3 <[Fe/H]< -0.5 



( Melendez et al.l (|2008f )) 



(iii) The thick disk low-metallic i ty ta il ends near 
[Fe/H]= -2.2 (jChiba and BeersI (|200(t l. 

(iv) Th e bulge forma t ion ti mescale is about 1 Gyr at 



most ( Zoccali et al 



(120081)). 

(v) The halo is made of two broadly overlapping struc- 
tural components: an inner halo with axis ratio 
close to 0.6, modest net prograde rotation, and peak 
metallicity at [Fe/II]= —1.6; and an outer halo 
with nearly spherical shape, net retrog rade rotation 
and p eak metallicity at [Fe/II]= —2.2 (jCaroUo et al 



(2003)). 

(vi) The mean rotational velocities are close to zero or 
negative for the outer halo {R ~ 12kpc), which may 
correspond to the radial h mit of the rapidly rotat - 
ing thick disk comp onent ( Chiba and Beers ( 2000h : 



ICarollo et al 



12003)). 

(vii) A linear dependence of mean rotational veloci- 
ties on the metallicity for [Fe/H]^ —1.6, appears to 
be consistent with the c ontraction of the inner halo 
( Chiba and Beers ( 2000l ) ) after an earlier assembling 
phase. Interestingly, the in ner halo metallicity peak s 



around [Fe/H ] ?« — 1 .6 (e.g..lRvan and NorrisI (|1991[ ): 



Carollo et aP (|2007h : IPrantzosI (| 200; 



(viii) The distribution of specific angular momentum 
appears to be consistent with both a halo-bulge 
and a thick disk-thin d isk collapse , but not with 



a spheroid-disk col l apse (Wvse and Gilmorel (1992) 
Ibata and Gilmorj (Il995l )). 



At this point, the question arises if the above men- 
tioned results are a natural consequence of the cur- 
rently favoured theory of galaxy formation based on hi- 
erarchical assembly of cold dark matter haloes in pres- 
ence of both baryonic matter an d quintessence (e.g.. 



iMota and van de Bruck ( 2004 ): Horellou and Bergg 
(|2005, )). In the framework of QCDM models, the initial 
density perturbations in the early universe have larger 
amplitudes on smaller scales, which makes low-mass 
density perturbations turn around and recollapse first. 



Accordingly, star formation is expected to start therein. 
For further details r efer to earlier a t tempts related to 
CDM scenarios (e.g., ICaimmil (|1997l ): IChiba and Beers 
(2003)). 

Low-mass density perturbations lying within a large- 
mass density perturbation, may be considered as sec- 
ondar y peaks bound to a primary peak (e.g., RvdenI 
( 19881 )). In particular, objects existing at present with 
Galactic mass are likely built up by the merging and 
accretion of smaller, less massive progenitors. It may 
safely be thought that these progenitors - at least in 
the denser regions of the universe - allow baryons to 
concentrate in their cores. This process takes place as 
an inevitable consequence of energy dissipation via ra- 
diative cooling and anelastic collisions, which does not 
carry significant amount of angular momentum. 

Stellar systems (and sufficiently clumpy gas clouds) 
moving along highly eccentric and/or highly inclinated 
orbits, make angular momentum transport possible via 
quadrupole interaction. More specifically, angular mo- 
mentum is transferred from the orbital motion of the 
tighly bound cores to the loosely bound outer halo par- 
ticles as the cores sink into the centre of the forming 
halo. This loss of angular momentum by the cores, due 
to dynamical friction, will lead to a low specific an- 
gular momentum in the assembled galaxy, as seen in 
ellipticals. F^or furt h er details refer to earlier attempts 
( Zurek et all (|l988l ): Ifiekki and Chibal (|2000l )). On the 
other hand, stellar systems formed from the gas heated 
in the course of merging and accretion, when a sub- 
stantial fraction of peculiar energy has been dissipated, 
would maintain their specific angular momentum, as 
seen in disk galaxies. 

In this view, the inner halo should have a fiat- 
tened density distribution with a fir iite prograde ro- 
tation , conformly to obse rvations (|Chiba and Beers 
( 200n^ : ICarollo et all (|2007h l and the angular momen- 
tum dist ribution may be similar to its bulge counter- 
part (Wvse and Gilmore ( 19921) ). The thick disk could 
have been ac creted from the extern a l part of the proto- 
Galaxy (e.g., Sommer-Larsen et al.l ( 1997t )). and its an- 
gular momentum distribution may be simila r to its thin 
disk counterpart (jlbata and Gilmorel (|l995l )). In other 
words, the inner halo and the bulge could take origin 
from the denser internal region of the initial density 
perturbation, where smaller spin growth via tidal in- 
teractions with neighbouring objects occurs, while the 
less dense external region, characterized by a larger spin 
growth, could be the progenitor of the thick and thin 
disk. The outer halo could be built up via accretion of 
smaller fragments or satellites even if, initially, bound 
to different peaks. 

The building blocks of each subsystem may be con- 
ceived in a twofold manner, as lower mass overdensi- 
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ties which first decouple from the universe (progenitors 
of compact dwarf galaxies) from the cosmological side, 
and upper mass gas media which form stars (molecular 
clouds or proto-globular clusters) from the astrophysi- 
cal side. 

In the framework of inhomogeneous models of chem- 
ical evolution considered in the current paper, star for- 
mation takes place within a fixed fraction of identical 
regions. At the end of any step, the enriched gas of 
active regions is homogeneously mixed with the unpro- 
cessed gas of quiescent regions, and a new set of identi- 
cal regions is produced, which makes the beginning of 
the subsequent step. Accordingly, the chemical enrich- 
ment history in different environments, is expected to 
be similar and even coeval for a strictly universal IMF. 
In fact, the total number of steps is an input parame- 
ter, but the duration of a step can freely be assigned 
(COO, COl, C07). For fiducial values of mass and radius 
related to the proto-bulge (M = Mio, R = 6kpc) and 
the proto-inner halo + disk (M = 6M10, R = 15kpc), 
the free-fall time is 0.08-0.12 Gyr, a few times lower 
than the b ulge formation time scale inferred from ob- 



servations ( Zoccali et al 



diooS)). 

Different environments could imply a different prob- 
ability of a region being active i.e. a different ratio of 
active to active -I- quiescent regions at any step and, in 
turn, different global efficiencies of the star formation 
rate. In this sense, chemical evolution could be influ- 
enced by dynamical evolution: for instance, a higher 
number of active regions corresponds to environments 
where tidal effects between neighbouring regions are 
strong (e.g., the bulge) and vice versa (e.g., the halo). 

Models implying large amount of gas inhibited from 
star formation (HI in the current paper) succeed in 
reproducing (incomplete) data related to the outer halo 
population. To this respect, some considerations seem 
useful provided the fit shall be confirmed using richer 
samples. Even if mass conservation formally holds in 
models under discussion, the chemical evolution acts 
as gas inhibited from star formation and a fraction of 
gas synthesized in short-lived stars and returned to the 
interstellar medium, were lost from the outer halo and 
presumably accreted by the proto-Galaxy. 

Using an outer halo star mass fraction, s/ = 0.13, in- 
hibited (i.e. lost) gas mass fraction, Df = 0.87, as hsted 
in Table O for model HI, and assuming a stellar outer 
halo upper mass limit equal to lO^m0, the total mass 
lost from the outer halo cannot exceed 6.7 lO^m©. The 
total mass of the proto-Galaxy is about ten times larger, 
which makes little influence on the related chemical evo- 
lution. In conclusion, mass loss from active regions 
makes a different oxygen enrichment with respect to 
environments (inner halo, thick disk, thin disk, bulge) 
where mass conservation holds to a better extent. 



The observed abundance distributions used in the 
current paper provide additional support to points (i), 
(ii), (iii), (v), listed above, where points (vi) and (vii) 
may be considered as complementary to (v) . The com- 
parison with model predictions, related to oxygen, is 
also consistent with the above mentioned points: the 
SN thick disk, the SN inner halo, and the bulge had 
similar and coeval oxygen enrichment within the com- 
mon abundance range. On the other hand, a different 
model with gas inhibited from star formation has to be 
used for the outer halo, which is consistent with points 
(v), (vi), and (vii). 

Keeping in mind point (iv), the above results imply 
a short formation timescale for both the (metal-rich) 
inner halo, the thick disk, the thin disk, and the bulge. 
Oxygen enrichment within a common range in different 
subsystems seems to be slightly dependent on dynam- 
ical parameters, according to point (viii). If any two 
subsystems formed on comparable timescales, as well 
as had a strictly universal IMF and star formation his- 
tories, the abundance pattern would be similar even if 
the two populations lack a physical connection. 

An idea of how galaxies did form is offered by celes- 
tial bodies which are presently assembling i.e. cluster 
of galaxies. Similarly, individual proto-galaxies may 
be conceived as made of subunits (Still visible dwarf 
spheroidals? Subsequently dead dwarf spheroidals with 
surviving nuclei in form of globular clusters? Globu- 
lar clusters?) which are virializing in the inner volume 
but are still falling in the outer volume, while a central 
structure (the proto-bulge) is growin g up. 



As outlined in an earlier attempt (|Chiba and Beers 
(l2000h ). more elaborate numerical modelling of the for- 
mation of large spiral galaxies, such as the Milky Way, 
is also needed in order to clarify the physical processes 
that lead to the currently observed dynamics and struc- 
ture of the spheroid and disk components. It is of par- 
ticular importance to model and understand the chemo- 
dynamical evolution of the system of sub-galactic frag- 
ments in the c ourse of the Galaxy assembling (e.g., 
Prantzosl(|l994l )). A fundamental understanding of the 
formation and the evolution of the Milky Way, is ex- 
pected to provide additional insights on the formation 
and the evolution of disk-type galaxies. 



6 Conclusion 

The main results of the present attempt may be sum- 
marized as follows. 

(1) The empirical differential oxygen abundance dis- 
tribution (EDOD) is determined from recent data 
involving different subsamples related to the solar 
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neighbourho od (SN) halo, thick d is k, thin disk, and 
to the bulge ( Ramirez et al. ( 2007t )): rMelendez et al 



imi)). 



(2) The EDOD related to different subsamples (includ- 
ing less recent ones) and/or different populations, 
show qualitative agreement above a threshold in iron 
abundance, [Fe/H]> —1, where the data in richer 
subsamples are statistically significant. The sole ex- 
ception concerns hal o stars, possibly due to subsam- 
ple incompleteness ( Ramirez et al. ( 20071)') an d /or 



different halo populations (| Car olio et al. 

(3) The SN thick + thin disk EDOD is deduced from 
its thick disk and thin disk counterparts, for differ- 
ent choices of the thick disk to thin disk mass ratio 
within the range, 0.1-0.9. The related effect is com- 
parable to the error calculated for the EDOD. 

(4) Inhomogeneous models with strictly universal ini- 
tial mass function (IMF), implying gas is neither 
inhibited from, nor enhanced in, forming stars at 
different rates with respect to a selected reference 
case, provide an acceptable fit to the EDMD related 
to different environments. To this aim, a SN thin 
disk low-metallicity tail (down to </> ~ 0.2) has to be 
assumed, and a single oxygen overabundant sample 
star must be disregarded as a peculiar object and/or 
affected by undetected systematic error. The exis- 
tence of a strictly universal IMF implies a similar 
chemical enrichment within active regions placed in 
different environments which, on the other hand, ex- 
hibit increasing probability of a region being active, 
and then an increasing global efficiency of the star 
formation rate, passing from the SN halo to the SN 
thin disk, the SN thick disk, and the bulge. The 
evolution of the SN thick -I- thin disk as a whole 



( HavwoodI (|2008[ )). computed using inhomogeneous 



models, is also consistent with the related EDMD. 
(5) A special case of stellar migration across the disk 
(inflowing stars with equal oxygen abundance with 
respect to the existing gas) can be described by in- 
homogeneous models with enhanced star formation. 
In particular, values within the observed range in 
SN gas mass fraction may allow a fraction of present 
SN stars were born in situ, and a comparable frac- 
tion is due to the net effect of stellar migration, ac- 
cording to recent results based on high-resolution A'^- 
body -|- smooth part icle hydrodynamics simulations 
iRoskar et all (|2008l )'l. 
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Appendix 
A Caption of symbols 

To help the reader, the meaning of the parameters hsted in the Tables throughout the paper, is explained below. 
The end of evolution is denoted by the index, /. For further details refer to parent papers (COO, COl, C07, COS). 

Zq Mass oxygen abundance. 

Mass oxygen abundance normalized to the solar value. 

A(j) Mass oxygen abundance (normalized to the solar value) bin, related to a selected [0/H] bin. 
(f)^ Upper (plus) and lower (minus) normalized oxygen abundance limit, related to an assigned normalized oxygen 
abundance bin, Acf). 

A^cj) Right (plus) and left (minus) normalized oxygen abundance interval semiamplitude related to an assigned 

normalized oxygen abundance bin, Acf). 
A*0^ Newly synthesised oxygen gas mass fraction (normalized to the solar value) within an active region, in 

inhomogeneous models of chemical evolution, at the end of a step. 
A*(p Newly synthesised oxygen gas mass fraction (normalized to the solar value) within the whole system, in 

inhomogeneous models of chemical evolution, at the end of a step during the expected evolution. 
ip Empirical differential oxygen abundance distribution (EDOD). 
Theoretical differential oxygen abundance distribution (TDOD). 
Upper (plus) and lower (minus) EDOD limit related to an assigned normalized oxygen abundance bin, A0. 
A^Tp Top (plus) and bottom (minus) EDOD interval semiamplitude related to an assigned normalized oxygen 

abundance bin, A</). 

Tpo Theoretical differential oxygen abundance distribution (exact value) at the starting point of the expected 
evolution (COl, C07). 

■(/'I (40) Theoretical differential oxygen abundance distribution (exact value) at the ending point of the first step 

of the expected evolution [COl, C07, Eq. (40)]. 
■(/'I (41) Theoretical differential oxygen abundance distribution (second-order approximation) at the ending point 

of the first step of the evolution in the general case [COl, C07, Eq. (41)]. 
/ip Gas mass fraction which allows star formation within an active region at the end of a step. 
D'^ Gas mass fraction which inhibits star formation within an active region at the end of a step. Negative values 

denote enhancement of star formation due to gas inflow. 
Long-lived star mass fraction within an active region at the end of a step. 
q Effective gas mass fraction within a region at the end of a step. 
^ Gas mass fraction which allows star formation. 

D Gas mass fraction which inhibits star formation. Negative values denote enhancement of star formation due to 

gas inflow, 
s Long-lived star mass fraction. 

a Mass fraction of a star generation which remains locked up in long-lived stars and stellar remnants. 
K Ratio of gas mass fraction which inhibits star formation to long-lived star and stellar remnant mass fraction. 
Negative values denote enhancement of star formation due to gas inflow. 
Lower mass limit of long-lived stars. 
p True oxygen yield. 

p' Effective oxygen yield related to both inhomogeneous star formation and inhibited or enhanced star formation. 
p" Effective oxygen yield related to inhibited or enhanced star formation. 
X Probability of a region being active. 



B Tabled] 

Table [8] is exceedingly large in the text, and for this reason it is broken to be completely accessible. 
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Table 12 Tabled right. 
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